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We have performed a global analysis of muon decay measurements to establish model-independent 
limits on the space-time structure of the muon decay matrix element. We find limits on the scalar, 
vector and tensor coupling of right- and left-handed muons to right- and left-handed electrons. The 
limits on those terms that involve the decay of right-handed muons to left-handed electrons are 
more restrictive than in previous global analyses, while the limits on the other non-standard model 
interactions are comparable. The value of the Michel parameter n found in the global analysis is 
— 0.0036±0.0069, slightly more precise than the value found in a more restrictive analysis of a recent 
measurement. This has implications for the Fermi coupling constant Gf- 
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PACS numbers: 13.35.Bv, 14.60.Ef, 12.60.Cn 



I. INTRODUCTION 

Muon decay, /j — > evv, is an excellent laboratory to in- 
vestigate the weak interaction. The energy and angular 
distributions of the electrons emitted in polarized muon 
decay are specified by the muon decay parameters p, rj, 
£, and 5 Q, 0] , conventionally referred to as the Michel 
parameters. The additional decay parameters £' and £" 
determine the longitudinal polarization of the outgoing 
electrons, and the parameters a, f3, a 1 , and (3' determine 
the transverse polarization Radiative muon decay, 
\x — » evuj, provides access to another decay parameter, 
fj Q. These 11 decay parameters, which are not all in- 
dependent, together with the muon lifetime provide a 
complete description of muon decay if the neutrinos are 
not observed. 

The decay parameters are related to the underlying 
space-time structure of muon decay by the most gen- 
eral, local, derivative-free, Lorentz-invariant transition 
matrix element, which can be written in terms of helicity- 
preserving amplitudes as Q|: 



v " 7 =S,V,T;e,iJ,=R,L 



^<e e |r> e ><^|r> M >. (i) 



non-zero in left-right symmetric models p|. For recent 
reviews of muon decay, see 0, 111 ■ 

The coupling constants g]^ represent 18 independent 
real parameters, so it is not possible to determine all of 
them from the existing muon decay observables. A global 
analysis of all the existing muon decay measurements 
can nonetheless set stringent model-independent exper- 
imental limits on the non-standard model contributions 
| llliA llllll l| that can then be compared to model calcula- 
tions such as those in |121j . which predict that g L ' R ' and 
9rl ' T snou ld be very small based on their contributions 
to neutrino masses. Recently, new measurements have 
been performed for the Michel parameters p and 5 
|l4j and for the transverse polarization parameters 77, 77", 
a', and f3' 15]. (rj and 77" parametrize the momentum- 
dependence of the CP-allowed transverse polarization in 
a manner different from, but equivalent to, a and f3. See 
Eq. JSJ) below.) Each of these new results is a factor of 
^2.5 more precise than the previous best measurement 
[l6| . This makes a new global analysis timely. In this 
paper, we present such an analysis to establish updated 
limits on possible non-standard model contributions to 
muon decay. 



In this equation, Gf is the Fermi coupling constant, de- 
termined from the muon lifetime, and the g]^ specify 
the scalar, vector and tensor couplings between /^-handed 
muons and e-handed electrons. Deviations from this ex- 
pression due to the non-local nature of the standard 
model weak interaction are 0(m^/m^) and, thus, are 
negligible compared to our current experimental knowl- 
edge of the muon decay parameters. Equation (|TJ is par- 
ticularly convenient because the standard model, with 
pure V — A coupling, implies g\ L = 1 and all the other 
coupling constants are zero. In contrast, additional cou- 
pling constants are non-zero in many extensions to the 
standard model. For example, g\ R , 9lr> an d 9rl are 
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II. FIT PROCEDURE 

Several different parametrizations for muon decay have 
been introduced over the years. One common version 
, which is based on the expression for the muon decay 
matrix element in charge-retention order, describes the 
muon decay observables in terms of 10 constants: a/A, 
a'/A, b/A, b'/A, c/A, c'/A, a/A, 0/A, a' /A, /3'/A. Each 
of these 10 constants is a bilinear combination of the 
coupling constants g]^ [8j. The normalization factor, A — 
a + 46 + 6c = 16, is fixed by adjusting Gf to reproduce 
the muon lifetime. The 9 remaining linearly independent 
terms may be determined from the experimental values 
of the muon decay parameters. This procedure was used 
in the global analysis performed in [!j ■ 

An alternative parametrization Q utilizes a different 
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set of bilinear combinations of the coupling constants: 



Qrr = 


1\9rr?1 


- \gU 2 , 


Qlr = 


\\9 S lr\ 2 + 


- \9Ir\ 2 + 


Qrl = 


\\9 S r L \ 2 + 


-\9 V rl\ 2 + 


Qll = 




\9 V LL \ 2 , 


Blr — 


^\9 S LR + 


§9 T lr\ 2 + 


Brl — 


Jq\9 S rl + 


§9 T rl\ 2 + 


la = 


1[9lr(9rl + ®9 RL 



_T |2 



_T |2 



= (a + ia')/2A, 



la = 



1[91l(9 S rrT + (9U*9L] 



= -2((3 + i{3')/A. (2) 

These bilinear combinations satisfy the constraints Q: 

< Qefj, < 1, where e, \i = R, L, 

0< B^ <Q efJ ., where e/i = RL, LR, 

\I a \ 2 < B LR B RL , 

\Ip\ 2 < QllQrr, (3) 



and the normalization condition, 

Qrr + Qlr + Qrl + Qll = 1- 



(4) 



The advantages of this parametrization are that it con- 
tains the maximum number of positive semi-definite bi- 
linear combinations of the gj^ and that the Q t ^ are 
directly interpretable as the total probabilities for /z- 
handed muons to decay into e-handed electrons 4] . Fur- 
thermore, experimentally it is found that Qll is close 
to unity. Thus, if Eq. (@J is used to eliminate Qll, the 
9 remaining variables are all close to zero. The bilinear 
combinations Qrr, Qlr, Qrl, B L r, B rt „ I „, h w ere 
adopted in the global analyses presented in jj, [T^J, IU • 

For convenience, in the present analysis we have 
adopted a hybrid set of independent variables: Qrr, 
Qlr, Qrl, B L r, B rl , a/ A, /3/A, a' /A, 0'/A. The 
muon decay parameters may be expressed in terms of 
these variables as 0, Hi : 



p = 


3 

i 4 


\(Qlr 


+ Qrl) - 


(Blr- 


VB RL ), 


s, = 


1 - 


2Qrr — 


w„ 

-^-Qlr + 


4 

2<9i?,L 


+ -^-(Blr — Brl 


£6 = 


3 

i ~ 


3 


7 

- ^Qlr + 


-^Qrl 


+ (B L r - Brl), 


= 


1 - 


2Qrr — 


2Qrl, 






e = 


1 - 


-^-{Qlr 


+ Qrl) ^ 


- -^{Blr. + B RL ), 



TABLE I: Experimental measurements included in the global 
analysis. 



Parameter 


Value 


Reference 


P 


0.7518 ± 0.0026 


\M 




0.75080 ± 0.00105° 


[13] 


5 


0.7486 ± 0.0038 


ii 




0.74964 ± 0.00130 


M 




1.0027 ± 0.0085 6 


\m 




0.99787 ± 0.00082* 


[12, 20] 


£' 


1.00 ± 0.04 


[16] 


e" 


0.65 ± 0.36 




V 


0.02 ± 0.08 


[16] 


a/ A 


0.015 ± 0.052 c 


[9] 


13/ A 


0.002 ± 0.018 c 


[9] 


'/ 


U.U 1 1 ± U.Uo 1 


|15| 


v" 


0.105 ± 0.052 d 


[15] 


a' /A 


-0.047 ± 0.052 e 






-0.0034 ± 0.0219 / 


[15] 


f3'/A 


0.017 ± 0.018 e 


m 




-0.0005 ± O.OOSO^' 


[15] 


Correlated with r) 


in ll '1. See text. 





Fit assumes P M = 1 . See text 
c p(a/A, /3/A) = -0.894 in [j]. 
d p\r),rf') = +0.946 in 
e p(a'/A, p'/A) = -0.894 in 
fp(a'/A,P'/A) = -0.893 in [1 



V = 

V = 

v" = 



-^(Qlr 



Qrl) 

(a - 2/3) /A, 
(3a + 2/3) /A. 



\{Blr 



Brl), 



(5) 



We have computed the joint probability density func- 
tion of the 9 independent variables using Monte Carlo 
integration techniques, in a manner similar to that de- 
scribed in For each experimental input, we have 
assumed that the corresponding probability distribution 
takes the form of a one- or two-dimensional Gaussian, 
truncated to the allowed parameter region if necessary. 
When statistical and systematic uncertainties have been 
quoted separately, we have added them in quadrature. 



III. INPUT VALUES 

Each of the previous muon decay global analyses 
0, 0, 0, E| utilized the same input parameters - p, 8, 
Pp^/P, a/ A, (3/A, a' /A, /3'/A. For each of these 

quantities, the previous analyses adopted the single most 
precise experimental measurement that was available at 
the time. 

We have utilized a different philosophy. We include 
all of the accepted muon decay parameter measurements 
that are reported in |T(| . with two exceptions. When 
|l6j | determines a decay constant from a single input, we 
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have included it. When |16| computes an average from a 
set of previous measurements, we have utilized this aver- 
age value in our fits. The exceptions involve the values 
adopted in [HI for r\ and the transverse polarization pa- 
rameters a/A, (3/ A, a' /A, and (3' /A. The adopted values 
for these five parameters are all derived from the global 
fit in [jj. Rather, for the transverse polarization parame- 
ters, we have used the experimental results |9J that were 
inputs to that previous global analysis. Meanwhile, we 
have not included the previous measurement of r\ in |l7j . 
Its uncertainty is a factor of 30 larger than the uncer- 
tainty in rj derived from our fit to the rest of the existing 
data. 

In addition to the measurements reported in |l6j| . we 
have included the recent measurements of the Michel pa- 
rameters p and 8 and of the transverse polariza- 
tion parameters T],T}", a /A, and (3' /A 0. The complete 
list of inputs is given in Table [I] 

Three of the input values require special consideration. 
In 0] > the value obtained for p assumed that r] is given 
by the accepted value, —0.007, and a contribution of 
± 0.00023 was assigned to the uncertainty in p associ- 
ated with the ± 0.013 uncertainty in the accepted value 
for j] 0]. We have performed two separate global fits. In 
one fit, we included the ±0.00023 in the uncertainty for 
p, then treated the measurement as uncorrelated with rj. 
In the other fit, we included the dependence of the mea- 
sured value of p on rj ^| an d reduced the uncertainty in 
p accordingly. The two fits give very similar results. The 
results for the latter are reported here. 

References 0] and |0, [2(J report measurements of 
and P^5/ p, respectively. In each case, P M repre- 
sents the polarization of the muon in pion decay. This is 
known to be > 0.99682 with 90% confidence Ha E3, E3l 
We have assumed P^ = 1 when including the mea- 
surement in our fits. In contrast, |l9ll20l | specifically men- 
tion the existence of possible unknown sources of muon 
depolarization prior to the decay and, therefore, choose 
to interpret the measurement of P^S/p as a lower limit. 
In recognition of this, we have performed two separate 
fits. In one fit, we assumed the probability distribution 
for £,5/p is Gaussian as given in Table |U which is equiv- 
alent to assuming P^ = 1. In the other, we treated the 
probability distribution as Gaussian when £,5/p is less 
than the central value and constant when £5/ p is greater 
than the central value. We report results from the for- 
mer fit because it finds slightly less restrictive limits for 
Qua, Qlr, and Blr, while the remaining parameters 
are unchanged between the two fits. 



IV. RESULTS 

Table ITT1 shows the results of the global fits. 90% con- 
fidence level upper limits (lower limit for Qll) are given 
for the positive semi-definite quantities. For complete- 
ness, mean and r.m.s. values are also specified for these 
variables, even though the corresponding probability dis- 



TABLE II: Results of the global fit. For the positive semi- 
definite quantitites, 90% confidence level upper limits (lower 
limit for Qll) are given, and for completeness, the mean and 
r.m.s. values from the fit are also given in parentheses. 



Parameter 


r ib r\,cbiiiL 1 a ±u 1 




<. 1.14 l^U.OU zt U.00J 


Qlr 


< l.y4 [l.ZZ zt U.OoJ 


Blr 


< 1.27 (0.72 ± 0.40) 


Qrl 


<44 (26 ± 13) 


Brl 


<10.9 (6.4 ± 3.3) 


Qll 


> 955 (973 ± 13) 


a/A 


0.3 ± 2.1 


P/A 


2.0 ± 3.1 


a'/A 


-0.1 ± 2.2 


(}'/A 


-0.8 ± 3.2 



tributions are far from Gaussian. In contrast, the output 
probability distributions for a/A, P/A, a'/A, and j3' /A 
are very close to Gaussian. 

Only a subset of the input quantities play a role in con- 
straining each of the fit parameters. Qrr is constrained 
primarily by the measurement of P^S/p. Qlr and B l r 
are constrained by the measurements of p, 5, and P^S/p. 
Qrl is constrained by the measurement of Brl is 
constrained by the combination of p, and 5. a/ A 
and a'/A are constrained primarily by the requirement 
I -^a 1 2 < B l rBr l . Finally, the large reduction in the 
uncertainties for a/ A and a'/A obtained when applying 
the I a constraint significantly reduces the correlations be- 
tween a/ A and j3/A and between a'/A and (3' /A. The fit 
finds p{a/A,(3/A) = -0.19 and p{a'/A, 0'/A) = -0.20. 
The correlations dominate the uncertainties in (3/ A and 
f3' /A quoted in Table [I] so the result is that the trans- 
verse polarization measurements become far more pre- 
cise. Overall, the constraints in Eq. © play a crucial 
role in reducing the uncertainties for a number of the 
fitted parameters, as was noted in 0. 

The fit results may be used to determine limits on 
the coupling constants g]^. These are given in Table 
IIIII Limits are also given for certain linear combina- 
tions of scalar and tensor interactions. Muon decay mea- 
surements alone are unable to separate contributions of 
g^ L from the standard model term, g\ L - Other experi- 
mental results, such as the branching ratio for 7r — > ev, 
argue that the scalar contribution must be very small, 
but this requires assumptions that go beyond the model- 
independent experimental limits that we are exploring 
here. In Q it was noted that the rate for inverse muon 
decay, i^e - — > fi~f e , may be used to distinguish between 
g£ L and g\ L - In Table ITTT1 the limits quoted for g^ L and 
9ll come from 8] , which utilized the inverse muon decay 
measurements in |23j. 

Table 1TTT1 shows that the present limits on |<?f o|, Iff^l, 
and |g^fj| are significantly better than those in y|, which 
arises from the inclusion of the new measurements of p 
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TABLE III: 90% confidence limits on the muon decay cou- 
pling constants in Eq. Q are compared to the previous ac- 
cepted values from |8j. Limits are also given for certain scalar- 



tensor interference combinations. 




Ref. [8J 


Present work 




< 0.066 


< 0.067 


\9rr\ 


< 0.033 


< 0.034 




< 0.125 


< 0.088 


\9lr\ 


< 0.060 


< 0.036 


\9lr\ 


< 0.036 


< 0.025 


\9rl\ 


< 0.424 


< 0.417 


\9rl\ 


< 0.110 


< 0.104 


\9rl\ 


< 0.122 


< 0.104 


\9 S ll\ 


< 0.550 


< 0.550 


\9ll\ 


> 0.960 


> 0.960 


\9lr + S9Ir\ 




< 0.143 


\9lr + *9Ir\ 




< 0.108 


\gl R - 2gI R \ 




< 0.070 


\9rl + 69rl\ 




< 0.418 


\9ltL + 2gl h \ 




< 0.417 


\9RL-2gl L \ 




< 0.418 



and 5. This was already noted in |14|. where calculations 
based on the measurements of p, 5, and P^S/p gave a 
slightly more restrictive limit for \gf R \ and a less restric- 
tive limit for |<7^j|- The new limits for and |<?J^| 
in Table ILTTl are slightly less restrictive than their previ- 
ously accepted values In @> these limits came from 
calculations in that were performed before an error 
in the reported value of P^S/p had been discovered. 
The present fit utilizes the corrected, less restrictive value 
of Pfj,t;5/ ' p |2ij. The new results for \gfi L \ and \g^ L \ are 
little different from the previously accepted values since 
the primary experimental input that constrains these, 
has not been remeasured since plj . 

The Michel parameter 77 merits special discussion be- 
cause it has an impact on Gp. Gf is typically calcu- 
lated from the measured muon lifetime assuming a pure 
V — A interaction. This leads to the quoted precision of 
AG F /G F = 9 x 10" 6 [3. However, a non-zero value of 
77 would change the muon decay phase space so that [2J| : 

G F *G%- A) (l-2rim e /m M ). (6) 



Thus, the ±0.013 uncertainty in the accepted value 
of r\ 0| leads to an additional uncertainty in Gp of 
AGf/Gf = 1.3 x 10~ 4 . The recent measurement of the 
transverse polarization of the electrons emitted in muon 
decay 15] included two analyses of the results. The gen- 
eral (model-independent) analysis led to the values of 77, 
77", a' I A, and j3'/A quoted in Table [I] A second, re- 
stricted analysis assumed that muon decay can be de- 
scribed with only two coupling constants, g\ L and g%pt- 
Equations J2J and (jSJ) imply that this assumption re- 
quires a/A — a'/A = and 77" = —77. The restricted 
analysis concluded 77 = -0.0021 ± 0.0070 ± 0.0010 and 
/3'/A = -0.0013 ± 0.0035 ± 0.0006. As noted above, 
the allowed ranges of a/ A and a'/A are severely con- 
strained by other muon decay data in the present global 
analysis. Thus, we find a model-independent result, 
rj = —0.0036 ± 0.0069, with slightly better precision than 
that of the model-dependent restricted analysis in |15| . 
This reduces the contribution of 77 to the uncertainty in 
G F to AG f /Gf = 6.7 x 10" 5 . Note that all of the new 
measurements [l3l Il4l 1 1 5| play important roles in reduc- 
ing the uncertainty in 77. 



V. CONCLUSION 

We have performed a new global analysis of all the ex- 
isting data on muon decay that do not include observa- 
tion of the outgoing neutrinos, including recent measure- 
ments of the Michel parameters p and 5 and the muon de- 
cay transverse polarization parameters 77, 77", a'/A, and 
f3' /A, The global analysis finds that the upper limits on 
the coupling constants Igf^J, Isl^I, an< ^ \9lr\ are niore 
restrictive than the previous accepted values. It also finds 
that 77 = — 0.0036±0.0069, which is nearly a factor of two 
more precise than the previous accepted value. 
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